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Abstract —Inter-cell interference (ICI) is one of the most 
critical factors affecting performance of cellular networks. In this 
paper, we investigate a joint link adaptation and user scheduling 
problem for multi-cell downlink employing HARQ techniques, 
where the ICI exists among cells. We first propose an approx¬ 
imation method on aggregated ICI for analyzing an effective 
signal-to-interference-and-noise ratio (SINR) with the HARQ 
technique at users, named identical path-loss approximation 
(IPLA). Based on the proposed IPLA, we propose a transmission 
rate selection algorithm maximizing an expected throughput at 
each user. We also propose a simple but effective cross-layer 
framework jointly combining transmission rate adaptation and 
user scheduling techniques, considering both HARQ and ICI. It 
is shown that statistical distribution of the effective SINR at users 
based on the IPLA agrees well with the empirical distribution, 
while the conventional Gaussian approximation (GA) does not 
work well in the case that dominant ICIs exist. Thus, IPLA 
enables base stations to choose more accurate transmission 
rates. Furthermore, the proposed IPLA-based cross-layer policy 
outperforms existing policies in terms of both system throughput 
and user fairness. 

Index Terms —Inter-cell interference, HARQ, link adaptation, 
user scheduling, cross-layer optimization. 

I. Introduction 

Compensation for uncertain fading phenomena is one of 
the most challenging issues in wireless communications. In 
order to improve link reliability and radio resource efficiency 
against such uncertainty, a hybrid automatic repeat request 
(HARQ) technique has been proposed in physical (PHY) 
layer m-m. Meanwhile, in medium access control (MAC) 
layer, dynamic link adaptation ii, a and user scheduling 
0, 0 techniques have been exploited by using channel 
state information (CSI) at the transmitter for point-to-point 
and multi-user environments, respectively. There have existed 
many studies on link adaptation considering HARQ for various 
fading channel models in point-to-point communications II- 
m . Moreover, several user scheduling algorithms considering 
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HARQ have been proposed for multi-user environments in 
single-cell networks 113 - 1321 . Through these studies consid¬ 
ering HARQ techniques, the resource efficiency of wireless 
communications has been improved. 

By the way, inter-cell interference (ICI) is another key 
factor to determine overall system performances in multi-cell 
networks. In traditional CDMA-based cellular networks in, 
the ICI was regarded as an additional source that deteriorates 
the performance in addition to intra-cell interference among 
users, which is typically managed by spectrum spreading 
(i.e., interference averaging) and power control techniques 
fl9l . In OFDM-based cellular networks like 3GPP LTE |20l . 
however, it has been observed that the ICI significantly de¬ 
grades the system performance, and thus, many techniques are 
being proposed in order to mitigate the ICI for OFDM-based 
cellular networks EH- Especially in heterogeneous network 
environments, there may exist dominant interferers which 
significantly affect adjacent cells l22l . Therefore, the ICI 
needs to be carefully managed in the OFDM-based multi-cell 
networks through efficient link adaptation and user scheduling 
algorithms. 

Recently, there have been several studies on HARQ-based 
multi-user systems in the presence of interference. Narasimhan 
analyzed throughput performance of the two-user interference 
channel with receiver cooperation | [23il . Denic proposed a 
robust HARQ-incremental redundancy (IR) scheme in the 
presence of unknown interference such as jamming IMS. 
For multiple-input multiple-output (MIMO)-based HARQ sys¬ 
tems, several HARQ techniques were proposed by taking 
into account inter-carrier interference and inter-antenna in¬ 
terference ED-IZZl- Racz et al. investigated an inter-cell 
interference coordination (ICIC) technique in the uplink 3GPP 
LTE system, considering HARQ techniques |28l . Makki et al. 
l29l proposed a coordinated HARQ scheme which reallocates 
the spectrum of a successfully transmitted user to a user 
requiring subsequent retransmissions in cooperative multi-cell 
networks. Shirani-Mehr et al. proposed an optimal scheduling 
algorithm based on game theory in a multi-user (MU)-MIMO 
system with HARQ technique in the presence of ICI |f30| . They 
investigated a joint optimization of user scheduling and trans¬ 
mit beamforming with HARQ in a distributed manner. To the 
best of our knowledge, however, there has been no such study 
that jointly investigates link adaptation and user scheduling 
with the HARQ technique in a multi-cell environment. 

In this paper, we investigate a joint link adaptation and user 
scheduling problem in a multi-cell downlink network, taking 
both HARQ and ICI into account. Main contributions of the 
paper are summarized as follows: 
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• A novel approximation model on the aggregated ICI 
at each user is proposed for enabling each base sta¬ 
tion (BS) to determine the optimal transmission rate, 
which assumes the dominant interfering terms at each 
user have identical path-loss statistics. Thus, we call it 
identical path-loss approximation (IPLA) method. The 
effectiveness of IPLA is examined by comparing it with 
the conventional approximation method, i.e., Gaussian 
approximation (GA). 

• An optimal rate selection algorithm with IPLA is pro¬ 
posed for maximizing the expected throughput of a single 
link in the multi-cell environment. Then, a simple but 
effective cross-layer framework is also proposed, which 
jointly combines link adaptation and user scheduling with 
the HARQ technique for the multi-cell environments. 

• The performance of the proposed cross-layer framework 
is evaluated in terms of cell throughput and user fairness 
through extensive system-level simulations. 

From the performance evaluation, it is shown that the well- 
known GA on ICI is not accurate in link adaptation and 
user scheduling with HARQ for multi-cell environments with 
some dominant interferers, while the proposed IPLA is highly 
accurate on the aggregated ICI and thus, it provides an efficient 
joint link adaptation and user scheduling policy. 

The rest of the paper is organized as follows. In Section HI1 
the system model is introduced. In Section [III] we propose an 
optimal link adaptation (transmission rate selection) algorithm 
for a single link by considering both the HARQ technique 
and the ICI. In Section IIV1 we propose a cross-layer frame¬ 
work jointly combining link adaptation and user scheduling 
and compares the proposed framework with the conventional 
strategies. In Section [V] we show the performance of the 
proposed framework in terms of cell throughput and fairness 
among users. Finally, we present concluding remarks in Sec¬ 
tion [yl] 


II. System Model 

Fig. HI illustrates the system model considered in this paper. 
We take into account a multi-cell downlink network where 
there exist (K + 1) base stations (BSs) with M transmit 
antennas and N users with a single receive antenna in each 
cell. Each BS is assumed to selects a single user for data 
transmission in this paper for simplicity. In Fig. E the BS 
in the center, called home cell, is denoted by superscript 
(0) and BSs in other cells are denoted by superscript (fc), 
fc £ {1,..., K}. Each BS selects a user within its coverage 
at each time slot (or scheduling interval), and transmits data 
with a random beamforming (RBF) technique which is also 
called opportunistic beamforming (OBF) ED, OH. As known 
in the literature, the RBF technique can achieve the system 
throughput with true beamforming when sufficiently large 
number of users exist in a cell, while it can significantly reduce 
the signalling overhead such as full CSI feedback for the true 
beamforming. 

The scheduled user in each cell receives a desired signal 
from its corresponding BS and the ICI signals from K other 
cell BSs. We focus on the user in the home cell without loss 



-> Desired signal channel 

- -> Interference channel 

Fig. 1. System Model 


of generality. The received signal of the scheduled user in the 
home cell (i.e., k = 0) is expressed as 
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y(o) = h (0) v (0) a; (0) + V 2 ;( fc )v (fc) x (fe) + n (0) (1) 

iJua ll u 0 V “o , »uo X U k ' n u 0 ! t 1 / 


fc=1 


where u k denotes the index of the selected user in the fc-th 
cell. hu„ and denote the channel vectors of the desired 
signal from the home cell (i.e., k = 0) and interference 
signal from the fc-th cell, respectively. v,, t . indicates the 
RBF vector for user u,k in the fc-th cell and Xy'J represents 
the source data symbol of user Uk in the fc-th cell, and 
rJuJ denotes the additive white Gaussian noise (AWGN), i.e., 
n-uj ~ CA/”(0, No) where No denotes the noise variance. 

In Eq. i[7]i. Iluq 1 and denote the MISO channel vectors 
including large-scale and small-scale fading components, i.e.. 
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(0) (k) 

where h u f m and gu 0 \ m denote the small-scale fading signal 
term of user u o from the m-th antenna of the home BS and the 
fc-th BS, respectively, and they are assumed to follow a circular 
symmetric complex Gaussian distribution with zero mean and 
unit variance, i.e., huo,m ~ CAfc(0,1) and g^ 0 ] m ~ CN(0, 1) 
for m = {1We assume a slowly varying channel 
condition and, thus, h u f m and gufm are quasi-static during a 
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single HARQ retransmission process. L^j and L\? 0 ’ denote the 
large-scale fading power terms regarded as path-loss of user uq 
from the home BS and the fc-th BS, respectively. is given 
by 10 • (do/duo^J , (k = 0,..., K) where PL 0 denotes 

the path-loss in dB at reference distance do, a denotes the 
path-loss exponent, and d^ represents the distance between 
user uo and the fc-th BS. When user and BS locations in 
the home cell are given by (r, 9) and (D( k \ip( k )), d^ 
can be calculated by T^^(fcrcos( 0 ^^Wy, 

( k = 0,..., K) where D = 0 (i.e., duj = r ). In general, 
since BSs are deployed at fixed locations in advance, the home 
BS can easily know the location information of neighboring 
BSs. Additionally, we assume that the home BS also knows 
the user location information through periodic measurement 
or feedback from the user. The 3GPP LTE system has been 
already supporting several user positioning methods even if 
global positioning system (GPS) signal is unavailable lf33l . 
Since we assume the perfect positioning in this paper, the 
performance can be degraded when the estimated user location 
is imperfect. The effect of position estimation error is beyond 
scope of this paper. 

Since the RBF scheme is considered as in 113 11 , the 
beamforming vectors in Eq. (Q3 are obtained by = 

(fc = 0 where [-] 7 de- 
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u k ,M 

notes the transpose of a vector, Vu k '. m = ffa, 
a m 6 [0,1], 9 m ~ Uniform]— 7 r, 7 r], and 
Ylm=i a,n = 1- Through a property of the RBF scheme, the 
second term in Eq. 0 , the sum of ICI terms, is derived by 
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where w$ ~ CJ\f(0, 1) and the last equality is derived from 
an isotropic property of complex Gaussian random variable 
0 . 

Consequently, the received SINR of the scheduled user in 
the home cell is 

s 

7 = 


x + i Ip 


(5) 


where s 

,k r (k) 

jk—1 


= || hi° 0 } vi? || 2 

„(fc )|2 


r (0) |„, (0) |2 

— -L/Uq \UJu 0 


X = 

J2k=i L ^ K o\ w u K o\ 2 , and p = where | 2 ] = P x for 

k = {0,..., K}. Here, is a known constant based on user 
location information at the home BS and we assume that the 
effective channel power gain of uo, | 2 , is perfectly known 
at the transmitter (i.e., home BS). It is reasonable because 
the home BS knows the RBF vector (i.e., v^,, ) for user uq 


in advance and we consider a quasi-static channel condition 
where it is possible to estimate the effective channel power 
gain perfectly for the desired signal channel (i.e., h!^). Hence, 
the desired signal power term s is a known constant at the 
transmitter. Furthermore, the inverse of transmit SNR term 
is negligible in interference-limited regime (i.e., high SNR 
regime). 

Throughout this paper, we consider the Chase combining 
based HARQ (HARQ-CC) protocol, in which every retrans¬ 
mitted information is same as the one at the initial transmis¬ 
sion. The HARQ-CC protocol is simple but obtains a sufficient 
benefit of HARQ from the combined power gain. Thus, it is 
widely used in practical wireless communication systems. 

III. Optimal Rate Selection for a Single Link 

In this section, we first mathematically formulate the 
effective SINR and delay-limited throughput (DLT) which 
represents an expected throughput under a given maximum 
allowable number of transmissions in HARQ-based systems 
ED, 0 , mi, (Ml- Then, we consider a well-known Gaussian 
approximation (GA) on the ICI with noise at users. Finally, 
we propose an identical path-loss approximation (IPLA) on 
the aggregated ICI term at users. We also obtain the optimal 
transmission rate maximizing the DLT for both approximation 
methods. 

A. Effective SINR and Delay-Limited Throughput 

First of all, the effective SINR after the n-th transmission 
attempt after HARQ-CC^ combining becomes 

n n 

7 <"> = £ 1i = £*TI 7 ? (6) 

2=1 

where 7 , represents the received SINR at the i-th transmission, 
Xi denotes the ICI power at the i-th transmission, which 
is expressed as Xi = XEi Luo \ w uo (U)\ 2 where fj indi¬ 
cates the time slot index of the i-th transmission, and s = 
Luo\wuo(U)\ 2 = T!°o ) |tt 4 o ) (£i )| 2 denotes the desired signal 
power. Since the RBF vector is kept during retransmissions 
and the desired signal channel is quasi-static during retrans¬ 
missions, s is a known constant for every (re)transmission of a 
single packet. In contrast, beamforming vectors (i.e., ) in 

other cells are independently varying according to scheduling 
decisions by other-cell BSs although interference channel 
vectors g „ 0 are quasi-static during their own retransmission 
processes. More specifically, new users can be scheduled 
after their own transmission successes in other cells during 
retransmissions in the home cell. This causes asynchronous 
scheduling among different cells, which implies that different 
cells suffer from different user scheduling instances. Accord¬ 
ingly, the aggregated ICI term is independently varying for 
every (re)transmission due to the independently varying other¬ 
cell beamforming vectors. 

Next, the distribution of the effective SINR needs to be 
analyzed for transmission rate selection. We start to derive the 

‘The HARQ-CC has been widely adopted in 3GPP HSPA (35), WiMAX 
(36), 3GPP LTE CD, and their evolutions. 
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distribution of the effective SINR based on numerical inversion 
of characteristic function from the following lemma. 

Lemma 1 (Inversion formula of Gil-Pelaez. 071? ): Let 
(j>(t ) = fy > e^ tx dF(x) be a characteristic function (CF) 
of the one-dimensional distribution function F(x). For x 
being the continuity point of the distribution, the following 
inversion formula holds true: 


1 1 

2 7r 


e~ itx <f(t) - ei tx (j)(-f) 


Xm 


2 jt 


dt 


dt. 


(7) 


where Im{-} denotes the imaginary part of a complex number. 
Proof: Refer to fl37l . ■ 


By using Lemma 1, if we know the CF of the effective 
SINR after the n-th (re)transmission which is denoted by 
0 7 ( n )(f), we can obtain the cumulative distribution func¬ 
tion (CDF), F 7 („)(x). Assuming the information-theoretic 
capacity-achieving channel coding scheme, the outage proba¬ 
bility after the n-th (re)transmission is defined by 

Pout(n, R) = Pr {log 2 (1 + 7 (n)) < R} = F j{n) (2 R - l) , 

( 8 ) 


where R denotes the required transmission source rate. Then, 
the DLT is obtained by lUTI 


N„ 


S(R)± Y -\Pout{i-l,R)-Pout(i,R) 

Z- J n 


i— 1 


l L 


(9) 


where N max denotes the maximum allowable number of trans¬ 
missions in an HARQ retransmission process. Substituting 
Eqns. 0 and 0 for Eq. 0, the DLT is finally rewritten 
by 


N max 


S(R ) = E 

2—1 


R 

i • 7r 



Xm\ 


0 -jt(2 R -l) 

t 


~ h(i- i)W) 

( 10 ) 


variable with shape parameter 1 and scale parameter s/a 2 z , % 
~ Inv-Gamma(l, s/cr|). 

First of all, the PDF of 7 % by the GA is 

fr i (x)=(s/a 2 z )x- 2 e- sl ^ x \ ( 12 ) 

To derive the distribution of sum of 7 \ (i.e., 7 (n)), the CF of 
7 i is derived first as follows: m 


07. (*) = 




(13) 


where K v (-) denotes the modified Bessel function of the 
second kind. Since 7 j’s are identically and independently dis¬ 
tributed (i.i.d.) random variables due to independently varying 
Xi, the CF of 7 (n) is obtained by 


11 <^ 7 i (t) 
2=1 




(14) 


By substituting Eq. d for Eq. d. the DLT of the conven¬ 
tional GA is obtained by 


Sga(-R) = 



i=i 
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i ■ it 



Xm^QAit] i, R) - 4' G A(f; i-l,R) 


dt, 

(15) 
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Finally, the optimal source rate for maximizing the DLT 
through the conventional GA is determined by 


R*g a = 


argmax 
R> 0 


Nrr 


1— 1 


R 


J Xm{'I>GA(t;i,R)-^GA(t-,i-l,R)} 


dt. 

(16) 


B. Link Adaptation with Gaussian Approximation (GA) 

Traditionally, the sum of ICI terms is widely approximated 
as a Gaussian distribution by the well-known central limit the¬ 
orem (CLT) for even six interference components considering 
7-cell structured cellular networks 081 . Il39l . Therefore, we 
investigate a rate selection scheme assuming that the ICI plus 
noise term follows a Gaussian distribution as a conventional 
link adaptation scheme. Through the GA, the effective SINR 
after HARQ-CC combining can be approximated by 


7 W = E 
2=1 


s 

Xi + 1 Ip 


Et< 


n 



(ii) 


where Zi denotes a complex Gaussian random variable with 
zero mean and variance of the sum of ICI and noise powers, 
i.e., Zi ~ CM ( 0 , er|) where <r| = YX-i Luo + 1 /p, and 
s = L^|ruio'*(fi)| 2 . Here, 7 , is an inverted Gamma random 


C. Link Adaptation with Identical Path-Loss Approximation 
(IP LA) 

Assuming the interference-limited regime (i.e., p 1), the 
effective SINR after HARQ-CC combining can be approxi¬ 
mated by 


7 (n) « E 

2=1 
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K 

k =1 


Lu k o\WuJ (U) | 2 


(17) 


where ti denotes the time slot index of the i-th transmission 
and s = Li^|wi^(fi)| 2 . The sum of ICI terms in the 
denominator of Eq. ( 1171 ) is a weighted sum of Gamma random 
variables since (ti) ~ CM(0, 1 ) and ^ Lul for 
k f l G {1 Note that there exist no closed-form 

expression for such distribution even though there have been 
some efforts to develop computational methods ED, El. 
Furthermore, the distribution of the effective SINR after the 
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n- th transmission attempt, 7(71), which is the sum of inverse 
of the weighted sum of Gamma random variables, has a 
much more complicated form and therefore it is intractable to 
derive its CF mathematically. Since only the contribution of 
the aggregated ICI rather than individual ICIs is interested in 
the effective SINR and even non-dominant interferers cannot 
be simply negligible^, we propose to approximate all path- 
loss terms from other-cell BSs to be identical as their average 
value. Then, the effective SINR can be approximated by 


7 (n) ~ = Y 1 

i—1 i= 1 


S 

E^Zlu^OI 2 ’ 


(18) 


where L = X ELi rf* 0 denotes the average value 
of all path-loss terms from other cells. It is worth 
noting that the proposed IPLA preserves the average 
statistics of the aggregated ICI since for a given user, 

E[Ef=iiko ) (*;)l 2 ] = ^ E [Efc=i \wuJ(U)\ 2 ] = 

£Ef=iE[l«& ) fe)l 2 ] = LK = Ef=i L$, While 

EfEti^oV^^)! 2 ] = T,k=i L uo E [\w^(u)\ 2 ] = 

Ef =1 L^uo because L^f’s are deterministic for the given 
user. 

Now, the sum of ICI terms becomes the sum of i.i.d. 
Gamma random variables and it also follows a Gamma 
distribution. After all, the approximated SINR at the i-th 
transmission, 7 *, follows an inverted Gamma distribution, i.e., 
7 i ~ Inv-Gamma(AT, s/L). Hence, the probability density 
function (PDF) of 7,; is 




( S / L ) K 

(K — 1 )! 


x~ K ~ 1 e~ s ^ Lx \ 


(19) 


Then, the CF of 7 * and 7 (n) are derived, respectively, as 
follows: l40ll 


fa, (t) = 


(K — 1 )! 


K k 


Ajst 
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{n){t) 11 07 i {t) 


2 (— 

(K — 1)! 


■K 
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( 20 ) 


( 21 ) 


where K v f) denotes the modified Bessel function of the 
second kind. By substituting Eq. (ED for Eq. d(, the DLT 
by the proposed IPLA is 


SlPLA(i?) = ]T 
2=1 


R 

i • 7 r 



Xm^ip LA (t;i, R) - ^i PLA (t;i - l,i?)} 


dt, 

( 22 ) 


where 


4'ipla (f; i, R) 


e -jt( 2 b -1) 


t 


H-fV 

{K — 1)! 


K k 



2 


2 To validate this statement, we examine the expected throughput of taking 
three dominant interferers among whole interferers, compared with the exact 
one in Section [V-BI 


Eventually, the optimal source rate based on the proposed 
IPLA for maximizing the DLT is determined by 


d* _ 

ii IPLA — 


argmax 
R> 0 



R 
i ■ 7r 



ipla (t; i , R) — ’Fipla^; * — 1, -R) j 


dt. 

(23) 


The above optimal source rate can be easily found by a grid 
search or a Golden section search with range between zero 
and a proper upper limit, since the DLT has a shape of quasi¬ 
concave function with respect to the source rate as shown 
numerically in Section [V-B I although it cannot be analytically 
proved due to a sophisticated form of the DLT formula. 


IV. Joint Link Adaptatin and User Scheduling: 
Cross-Layer Framework 

A. Overall Procedure 


We first propose a simple cross-layer framework to per¬ 
form both link adaptation and user scheduling considering 
HARQ and ICI. The cross-layer framework consists of three 
main components: rate selection , effective rate mapping, and 
scheduler. The roles of components are described as in the 
following. 

1) Rate Selection (RS): The RS plays a role to determine 

an optimal transmission source rate for the u-th user 

at initial transmission instance of the HARQ-based system. In 
this paper, we consider RS schemes to maximize the DLT of 
each user considering HARQ retransmission and ICI statistics 
as presented in the previous section. 

2) Effective Rate Mapping (ERM): The transmission source 

rate is different from the achievable rate in HARQ-based 
systems due to uncertain retransmissions. Therefore, an ef¬ 
fective rate, which is close to the achievable rate, needs to 
be taken into account for user scheduling if it is available. 
The ERM determines an effective rate Rctr.u(t) for the u- 
th user as a function of the optimal source rate i.e., 

Re&,u(t) = f(R^{t)), in order to adjust the scheduling priority 
of each user. Through such ERM, the instantaneous rate R u (t) 
in the scheduler is replaced by the effective rate /Crr.uff). After 
all, the scheduler selects a user with the highest utility value 
substituted into the effective rate. 

3) Scheduler: The scheduler determines which user is the 
best at every scheduling instance. There are three representa¬ 
tive scheduling algorithms: Round Robin (RR), Max C/I, and 
Proportional Fair (PF). In this paper, we take into account 
the PF scheduler for an asymmetric user distribution scenario 
where users have different distances from the BS, in order to 
consider user fairness. The PF scheduler is simply expressed 
as: 


* 

u 


= argmax 

■wElI 


Ru{t) 
Tu(t) ’ 


(24) 


where n denotes the set of users in a cell, R u (t) denotes 
the achievable rate of the w-th user at time slot t, and T u (t ) 
denotes the average throughput of the u-th user at time slot t, 
which is updated as T u (t + 1) = (1 — p) • T u (t) + j- ■ R u (t) ■ 
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I{« = u*} where t, denotes the pre-determined windowing 
interval for moving averaging and !l{x} denotes the indication 
function which is one if x is true and zero otherwise. In the 
HARQ-based systems, R u (t) should be modified considering 
the HARQ retransmission process and it can be done by the 
ERM in this cross-layer framework. 

The operating procedure according to the proposed cross¬ 
layer framework is illustrated as follows: 

(Step 1) [Rate Selection]: Determine R* u (t) 

(Step 2) [Effective Rate Mapping]: 

Determine Res >u {t) = / (i?*(t)) 

(Step 3) [Scheduler]: Determine u* = argmax 

(Step 4) [HARQ Transmission] 

- u* transmits with until successful transmission 

or maximum transmission limit. 

- Go to (Step 1) for all users after the end of the 
(re)transmissions of the scheduled user u*. 


B. Proposed Cross-Layer Policy and Other Candidates 

In this subsection, we propose an IPLA-based cross-layer 
policy. We also introduce two reference and two conventional 
policies for the performance comparison in next section. 
Hereafter, each cross-layer policy is denoted by the rate 
selection and effective rate mapping, i.e., V{ RS, ERM}, since 
all policies employ the same PF scheduler. 

1) Proposed IPLA-based Policy, V {RS-IPLA, 
£ipla(-Ripla)} ; The proposed IPLA-based policy is 
based on RS through the IPLA on the aggregated ICI term 
(so called RS-IPLA). According to the RS-IPLA proposed 
in the previous section, the transmission source rate of the 
w-th user is determined by Eq. dH. Since it takes into 
account an HARQ retransmission process using statistics of 
the aggregated ICI term, which is assumed by a Gamma 
distribution through the IPLA, the achievable rate is different 
from the transmission source rate. Thus, the expected 
throughput -based ERM is considered for user scheduling as 


Reff,u, — 'S'lPLA(^?IPLA,ix) 

N„ 


= E 


R 


IPLA,u 


i=l 


Jm|4>IPLA(f ; i, -ft IPLA,ix) 


^ $IPLA(i;*- 1,-RlPLA.ix)} 


dt, (25) 


where 


'TiPLAO; L -RlPLA.il) 

e -jt{ 2 *fPLA,u_ 1 ) 

t 


*Hr)‘ 

(*-1)1 


L if 



Based on the selected transmission source rate and effective 
rate, the user scheduling is performed with the PF criterion and 
then, the HARQ transmission is performed for the scheduled 
user according to (Step 4) in the operating procedure of the 
cross-layer framework. 


2) Reference Policies: 

• Genie-Aided Policy, 'PjRS-Opt, R ( )pt } 

The genie-aided policy has perfect knowledge on in¬ 
stantaneous ICI terms. In this case, the transmitter can 
accurately adapt to instantaneous interference channel 
and the channel capacity depending on the instantaneous 
interference channel conditions is achieved without any 
retransmission and outage. Even though this policy is 
rather unrealistic, it offers an upper bound of the system 
performance. According to the RS-Opt scheme with the 
perfect knowledge for interference channels, the trans¬ 
mission source rate of the u -th user is determined by 


R *0 Pt,u 


log 2 1 + 


Eti L 


( fc )i„„( fe ) 


w^’ihW + l/p, 


(26) 

where t\ denotes the time index at initial transmission, 
|wi fc ' ) (ti)| 2 represents the exact effective interference 
channel power gain from the k -th BS to the n-th user 
in the home cell, s = L^ |wi°' ) (ii)| 2 , and p = 
Next, since the genie-aided policy does not cause outage 
and retransmission, the instantaneous rate-based ERM is 
considered as 


Reff.u — R, 


Opt, if' 


(27) 


• Instantaneous SINR Policy, 'P{ RS-/-SINR, i?* gINR } 

The instantaneous SINR policy is the simplest one based 
on the SINR value fed back from the receiver. This policy 
has inaccurate rate selection due to independently varying 
interference channels for every (re)transmission and a 
feedback delay. According to the RS-/-SINR scheme 
using the outdated feedback channel information, the 
transmission source rate of the vi-th user is determined 
by 


-Ri-SINR,u — 1°§2 


1 + 


Ef=i4 fe) k fc) (G-<5)| 2 + l/p, 


(28) 

where |u4 (Q — <5)| 2 represents the interference power 
gain from the //-th BS to user u with the feedback delay S, 
s = Li°' , |wi°' ) (ti)| 2 , and p = Since the instantaneous 
SINR policy takes advantage of only instantaneous infor¬ 
mation without consideration of HARQ retransmission, 
the instantaneous rate-based ERM is also considered as 


Reft, u 


= R 


i-SINR.ix- 


(29) 


3) Conventional Policies: 

• Average Interference Policy, 'P{RS-Avg-A’, } 

The average interference policy exploits an average value 
for the aggregated ICI term (26), since each ICI term is an 
uncertain and independently varying factor. According to 
the RS-Avg-A scheme replacing the aggregated ICI term 
by the average value, the transmission source rate of the 
u-th user is determined by 


R 


Avg -X ,u 


log 2 



(30) 


X + l/p 
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where X = E 


E 


K 

k= 1 


4fW } | 2 


E 


fc=l 


L 


(*) 


s = 


and p = The average interference 
policy has an identical source rate during retransmissions 
because of applying the average value for the sum of 
ICI terms. Thus, similarly to Eqns. (l27l> and (l29l >. the 
instantaneous rafe-based ERM is considered as 


Heff ,u — ^Avg -X,u' (31) 

• GA-based Policy, T^RS-GA, 5'ga(^?ga)} 

The GA-based policy is based on RS through GA for 
the sum of ICI and noise terms. According to the RS- 
GA scheme investigated in the previous section, the 
transmission source rate of the u-th user is determined by 
Eq. (IT6t . Since the GA-based policy considers an HARQ 
retransmission process using statistics of the aggregated 
ICI term, the expected throughput- based ERM is consid¬ 
ered as 


N maa 

R e ff,u = S G a(Rga,u) = ^2 

i= 1 


E>* 

n GA,u 
i • 7r 



Xm|$GA(i; i, R*ga,u) - $ga( t; i 


1 ,R, 


'GA,u) } 


dt, 

(32) 



Fig. 2. Q-Q plots of the proposed IPLA and the conventional GA (K = 6, 
r = 250 m, 9 = n/2, a = 3, |u!^ | 2 = 1, Nm ax = 4) 


interference-limited situation, we set transmit SNR, p, to 43 

dB0 


where 


$GA(i;*,-RGA,J 



V. Numerical Results 

In this section, we first examine the effectiveness of the 
proposed IPLA through a quantile versus quantile (Q-Q) plot 
on the effective SINR distribution, F^t n \(x). After that, we 
discuss effects of user distance and path-loss exponent on 
the rate selection in a single link. Finally, we evaluate the 
performance of the proposed, conventional, and reference 
cross-layer policies in terms of system throughput and fairness 
metric, through system-level simulations. As basic simulation 
setups, we consider an 1-tier cellular network with six other 
cells (i.e. K = 6) where users are asymmetrically distributed 
in the home cell. The BS-to-BS distance is set to 1000 m (i.e., 
£)( fc ) = 1000, Vfc) and angles between BS in the home cell 
and BSs in the other cells are set to ^ i.e., 

■if = §i f > — f j — f j — ^p]. We set the distances between 

users and the home BS to r £ [150, 200, 250,300,400]m and 
each element in the vector is equally set according to the 
number of users. Therefore, we just consider the number of 
users as multiple of five and it is 250 m for a single user case. 
Additionally, the angles between users and the home BS are 
uniformly determined as 6 =Uniform[— n, 7r]. For path-loss, 
we set PL 0 to 37 dB at reference distance do = 1000 m and 
path-loss exponent, a, to 3. The maximum allowable number 
of transmissions, N max , is set to 4, which is a typical value 
in LTE and WiMAX systems. In order to take into account an 


A. Statistical Distribution of Effective SINR based on IPLA 

In order to examine the effectiveness of the proposed IPLA, 
we introduce a Q-Q plot which is widely used for quanti¬ 
tative comparison between two distributions. It can provide 
an intuitive comparison between two statistical data sets as 
well as two theoretical distributions and more information 
on the local agreement between two distributions than other 
fitting tests such as Chi-square and Kolmogorov-Smirnov tests 
l43l . In this paper, we compare two theoretically approximated 
distributions with the real empirical distribution. In the Q-Q 
plot, the x-axis is based on the theoretical distribution with the 
approximated CDF, which is obtained by inverting the CDF, 
F~/ n ){ x )> an d the y-axis is based on the empirical quantile 
from a sample data set on the effective SINR obtained by 
statistical realizations. 

Fig. [2] shows the Q-Q plots of the proposed IPLA and the 
conventional GA compared with the real empirical distribution 
according to the number of transmission attempts of a single 
packet, n. Since the line y = x represents the identity of 
two compared distributions in the Q-Q plot, the proposed 
IPLA almost following the y = x line agrees well with 
the real empirical distribution, regardless of n values. As a 
general trend, with increasing the number of transmission 
attempts, n, the Q-Q plots of both approximated distributions 
move to the right-upper side, which implies a larger effective 
SINR value due to an HARQ-CC combining gain. While 
the Q-Q plot of the proposed IPLA agrees well, that of the 
conventional GA is flatter than the line y = x. This implies 
that the approximated distribution by the GA is more dispersed 
than the real empirical distribution. Additionally, difference 

3 In this setting, the average received SNR without interference becomes 6 
dB when the distance is 1000 m. 
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between the approximated distribution by the GA and the 
empirical distribution increases as the value of effective SINR 
increases. Through the comparison of the approximated distri¬ 
butions with the real empirical distribution, it is shown that the 
proposed IPLA offers a good approximation on the effective 
SINR, while the conventional GA gives significant differences 
in the approximation. 


B. Link Adaptation for a Single Link: Effects of User Distance 
and Path-Loss Exponent 

In this subsection, we investigate the effects of user distance 
and path-loss exponent on the rate selection according to the 
proposed and conventional link adaptation schemes. Fig. [3 
shows the DLT for varying source rate It in three different 
user distance values. Basically, the DLT has a shape of quasi¬ 
concave function and a single optimal point with respect to 
the source rate. As the distance decreases, a higher DLT is 
achieved since the desired signal power increases while a 
closer distance to the BS fundamentally yields smaller interfer¬ 
ence from other-cell BSs. In the figure, the solid lines denote 
the exact simulation results with perfectly known individual 
ICIs. Additionally, the dotted lines denote the simulation re¬ 
sults with perfectly known three dominant ICIs, which neglect 
the other three ICIs. Compared to both simulation results, 
neglecting non-dominant ICIs yields over-estimated DLTs due 
to the reduced interference even if it shows similar shapes 
of curves. On the contrary, the DLT analytically derived by 
the proposed IPLA has a high similarity with one by the 
exact simulation for all distance values. After all, the optimal 
source rate determined by the proposed IPLA is approximately 
identical to the actual optimal source rate on the exact simula¬ 
tion curves, regardless of the user distances. However, optimal 
source rates based on the two conventional link adaptation 
schemes, GA and average interference schemes, exhibit sig¬ 
nificant differences from the actual optimal value on the exact 
simulation curves. The gap between the optimal source rate by 
the proposed IPLA and one by the conventional GA increases 
as the user distance increases (i.e., as ICI increases), whereas 
the gap between the optimal source rate by the proposed 
IPLA and one by the conventional average interference scheme 
increases as the user distance decreases. Therefore, we can 
conclude that the conventional GA is relatively good for near- 
BS users while the conventional average interference scheme 
is good for edge users, in the link adaptation perspective. 

Fig. [4] shows the DLT for varying source rate R in three 
different path-loss exponent values. As the path-loss exponent 
value increases, a higher DLT is achieved since the interfer¬ 
ence is reduced with increasing the path-loss exponent value. 
The optimal source rate by the proposed IPLA also agrees well 
with the actual optimal value on the exact simulation curve, 
while those by the conventional schemes show significant 
differences. The basic trends of the differences are similar 
to those in Fig. [3] Consequently, for the conventional GA, 
the more interference exists, the larger difference occurs in 
the optimal source rate, while for the conventional average 
interference scheme, the less interference exists, the larger 
difference occurs in the optimal source rate. Fundamentally, 



Fig. 3. Effect of User Distance on Rate Selection (K = 6, 9 = it/ 2, a = 3, 
|t4° 0 ) | 2 = 1. N max = 4) 



0123456789 

Source Rate, R 


Fig. 4. Effect of Path-Loss Exponent on Rate Adaptation (K = 6, r = 250 
m, 6 = 7r/2, \Wuq\ 2 = 1 , Nmax = 4) 


both the conventional schemes exhibit significant differences 
with respect to the optimal source rate for large SINR values 
(i.e., a = 4), which correspond to the weak interference 
situation. However, in a scheduling-based multi-user system, 
a user with a large SINR value has more opportunities to 
be selected as the best user. Therefore, it is expected that 
the proposed IPLA-based cross-layer policy can obtain a 
significant throughput gain in the viewpoint of both link 
adaptation and user scheduling, compared to the conventional 
policies. 


C. System-Level Performance Evaluation: Cell Throughput 
and Fairness 

Fig- 0 shows the system performance of various cross-layer 
policies for varying the number of users in the home-cell. 
Specifically, Fig. 0 (a) shows the system throughput of the 
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proposed, conventional, and reference policies. The genie- 
aided policy, P{RS-Opt, f?Q pt }, provides an upper bound of 
the system throughput even if it is unrealistic. The instan¬ 
taneous SINR policy, ^{RS-i-SINR, /f* S | NR }, achieves the 
worst system throughput due to rather inaccurate estimation 
of the ICI term caused by the channel feedback delay. The 
proposed IPLA-based policy, ^{RS-IPLA, 5 'ipla(^?ipla*)}, 
always outperforms the conventional policies, V { RS-GA, 
Sga{Rqa)} anc * T’jRS-Avg-A'’, R%, & . x }> i n the entire range 
of the number of users, while both the GA and average inter¬ 
ference policies achieve almost identical system throughput. 
Note that although the GA-based policy exploits statistics 
of the ICI term, it achieves similar system throughput to 
that of the average interference policy which just utilizes an 
average value of the ICI term. Moreover, it achieves rather 
smaller system throughput than that of the proposed IPLA- 
based policy which also exploits equivalent average statistics 
of the ICI term. 

Fig- 0(b) shows the fairness metric performance of various 
cross-layer policies for varying the number of users in the 
home-cell. We consider the fairness metric in 03, G2, GO, 
1441 defined as TM (T),..., T N ) = l°g(T u ) where 

T u denotes the achieved throughput of the u-th user and 
N is the number of users in the system. As investigated in 
the previous work, the fairness metric offers a performance 
measure considering both system throughput and user fairness 
together. Under the PF scheduling algorithm with averaging 
time scale t c = oo, the fairness metric is maximized almost 
surely among the class of all schedulers ED- As shown in 
Fig. [5] (b), the proposed IPLA-based policy, V{ RS-IPLA, 
£ipla(-Ripla)}’ a ^ so outperforms the other three policies 
except for the genie-aided policy for all the number of users. 
The average interference policy, 'P{RS-Avg-A’, R* nst }, rather 
outperforms the GA-based policy, V { RS-GA, Sc a ( Rq A )}, 
in terms of the fairness metric, even if the GA-based policy 
exploits more information for ICI term than the average 
interference policy. It comes from the fact that the GA on the 
aggregated ICI is inaccurate when there exist some dominant 
ICIs which are general in OFDM-based cellular networks, 
although the GA on the aggregated ICI is well-approximated 
for a sufficiently large number of independent and identically 
distributed interferes. In contrast, the proposed IPLA is highly 
accurate in this environment. Accordingly, the proposed IPLA- 
based cross-layer policy is able to be the most efficient in 
the OFDM-based cellular networks where there exist some 
dominant ICIs. 

VI. Conclusion 

In this paper, we investigated a joint link adaptation and 
user scheduling in multi-user and multi-cell environments, 
considering HARQ techniques. Based on the proposed mathe¬ 
matical approximation method for the ICI signals the optimal 
transmission rate selection algorithm in terms of the expected 
throughput is proposed. As for multi-user environments, a 
novel and effective cross-layer framework combining the link 
adaptation and user scheduling is also proposed. Through 
extensive link-/system-level simulations, it is shown that the 



Fig. 5. System Performance in an Asymmetric User Distribution Scenario (a) 
System DLT vs. Number of Users (b) Fairness Metric vs. Number of Users 
(K = 6, r € [150, 200, 250, 300, 400]m, 9 =Uniform[-7r, vr], a = 3, 
Nmax — 4 ) 

proposed cross-layer policy significantly outperforms the con¬ 
ventional policies in terms of both cell throughput and user 
fairness. With consideration of both HARQ and ICI, we have 
tried to investigate more general and practical communication 
scenarios including multi-user MIMO, receiver beamforming 
at users with multiple receive antennas, and HARQ technique 
with incremental redundancy, but they were not mathemati¬ 
cally tractable unfortunately. Therefore, we leave these issues 
for future work. 
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